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Abstract 

We calculate, in the Glauber-Mueller approach, the ratio of the diffrac- 
tive dissociation cross section to the total cross section on nuclei. We ob- 
serve that shadowing corrections ( SC ) are significant, in the calculated 
a'^'^^^ jotot ratio but they do not lead to a smooth energy dependence of 
the ratio for a proton target as was observed at HERA. 
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1 Introduction 



In this letter we discuss the energy dependence of a'^^-^-^/crtot, where cr'^*-^-^ is the 
cross section for the diffractive dissociation of a virtual photon (7*), and atot is 
the virtual photon nuclear total cross section. It was observed at HERA that 
for DIS on a nucleon this ratio, when taken for definite masses of diffractively 
produced hadrons, is almost independent of W (the cm. energy of the vir- 
tual photon - proton ) , over a wide range of photon virtualities H]. At 
present there is no theoretical explanation of this intriguing phenomena, except 
for the quasi-classical gluon field approach 0. However, in this model the to- 
tal and diffractive cross sections do not depend on energy, in contradiction to 
the experimental data. Some attempts have been made to develop a successful 
phenomenology, which describes the experimental data fairly well Recently, 
aiming to explain this phenomenon theoretically , Yu. Kovchegov and L. McLer- 
ran have suggested that shadowing corrections (SC) play a crucial role in 
the kinematical region of interest. They derived a formula for a'^^-^-^ /atot which 
takes into account the quark - antiquark final state. In the following, we gen- 
eralize the Kovchegov-McLerran formula using the Glauber-Mueller approach 
to SC. We take into account the qq plus one extra gluon final state without 
which we are unable to explain the diffractive production at HERA p. In the 
same spirit we derive a formula for cr'^^^ ^ / cr tot for an arbitrary nucleus. 

For the case of a nucleus target, the Glauber-Mueller (G-M) approach can be 
justified theoretically P, P], while for DIS on a nucleon this approach is just one 
of the many formalisms used to estimate shadowing corrections. The kinematic 
region where the G-M approach is valid is well defined. In addition, as shown by 
Kovchegov , it provides the initial conditions for the GLR evolution equation 
[10[] in the full kinematic region of DIS. Q 

In section 2 we calculate cr'^*-^-^/ cxtot for a nucleus with qq and qqG final states. 
In section 3 we give numerical estimates for different nuclei and different virtual- 
ities Q^. In section 4 we conclude with a summary where we list some possible 
experimental consequences of our calculation. An early rough calculation of 
diffraction off nuclei has been published in Ref.p]. 



^ The theoretical proof assumed sufficiently heavy nuclei. Based on much experience with 
"soft" Pomeron models [O , it appears that the Glauber approach can be used for nuclei with 
A > 30. 
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2 Shadowing Corrections 

It was shown in that the total cross section for the interaction of a virtual 
photon with the target can be written as 

o-totil*+p) = J dz J (fr±P'^*{z,r±;Q^)adipoieixB,r_L) (1) 

= 2 J jdzj dWP^\z,rr,Q')lma'',\^,i^{xB,rr,b^), 

where P^* {z, r±; Q"^) denotes the probability of finding a quark-antiquark pair of 
size r_|_ inside a virtual photon, and (Tdipoie is the cross section for the interaction 
of a colour dipole with the target. The optical theorem has been used to express 
atot in terms of the forward elastic amplitude a%poif,- The s-channel unitarity 
constraint reads 

2Ima^'(s,6x) = + ^^"(s, 6^), (2) 

where G*" corresponds to the contribution of all inelastic channels. The total, 
elastic and inelastic cross sections of the process are given by 



atot = 2j d%i_lma'\s,b^) (3) 
aei = j d%^\a-\sM\'' (4) 

CT^n = j d%^G''^{s,h^). (5) 

Using these relations, the cross section for diffractive production of a qq pair is 

.^^ff^^^p^qq) = J d'b^P^\z,r^-Q')K\^^,^{xB,rr,b^)\^ (6) 



a 



To obtain an estimate for the amplitudes a*^' and G*" we assume that the 
dipole-target amplitude is purely imaginary, this is a reasonable approximation 
at high energies. Thus, the general solution of the unitarity constraint, given 
by Eq.(ED can be written 



a^'(x,rx;6x) = ^ [l - e (7) 

G7^poie{^,rr,b^) = i-e-^''(-'^-M_ (8) 

In the spirit of the Glauber-Mueller ^ approach, one may consider the 
function as an amplitude for the exchange of a hard Pomeron (gluon ladder), 
which is given by ( see Fig. 1) 



where S{b±) is the nucleus profile function for which we take a Gaussian form 
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Sib±) = ^e (10) 



In the Woods-Saxon parametrization |15], the size of the nucleus is given by 



Ra{WS) = 1.3 fm A^/^, while in the Gaussian parameterization R\ = 2/5R\{WS). 
A full explanation of all factors in Eq.(lO) are given in Refs. |^ and |T^. 

In the Glauber-Mueller approach at large Q^, a useful approximation is to 
replace the integration over z as follows 0, 



/O V ' -'-^ ^ 3^2Q2 

which reproduces the Kovchegov-McLerran formula 0, |T3| 



;i2) 



There are two kinds of corrections to this simple picture [|T3| (see Fig.l): 
1) As has been discussed 0, the G-M formula ( Eq.( |^) describes the rescattering 
of the fastest dipole interacting with nucleus ( see Fig.l-a ). This rescattering 
depends on the opacity fl^, given by Eq.(9). However, Eq.(lO) depends on 
AxG^*^^"^^ which is also screened due to the rescatterings of the fastest gluon 
in the gluon ladder (see Fig. 1-b ). To calculate this correction we iterate 
Eq.( 1^), substituting instead of fl^ the expression from the Glauber-Mueller 
formula ^ for rescattering of a gluon in a nucleus. 

= l7r'asArlxG^^^^P{x,^)S{b^) (13) 
^n^'^ = \7c'asrlxG^''ix,4^;b^) , 



where 



(^>rr;M = -7/ -7 L T7f2h-e ^ • (14) 



2 ) 9 / I I ^ '4 

rj_ TT^ Jx X Jr]_ rj^ y y 

Note that the gluon ladder gets an additional colour factor VLq = |^^- 

So far we have assumed that the diffraction dissociation final state is merely a 

quark-antiquark pair plus the original nucleus { qq + A state ). 

2) The second improvement is a correction to this final state. For example, 
this can be an excited nucleus or the emission of an extra gluon (see the cuts in 
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Fig.l-b 



Figure 1: Total cross section and diffractive cross sections for a dipole - nucleon 
interaction in the Glauber- Mueller approach (Fig. 1-a) and in the first inter- 
action of the G-M approach. The dashed lines show the diffraction dissociation 
cuts. 
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Fig.l-b ). In Ref. we incorporate these correction using the AGK cutting 
rules [Q. The inelastic cross section of a dipole scattering off the target with 
the emission of one extra gluon in the final state is given by 



^Zoieim - qqG) = I d%^e~^''"^^---Ma'^i( , (15) 

where 

diffir, ^ Nc^sA f^dx' [dr'l( -na^^^^y^V 
-oiIib,) = —^l-j^[l-e . j, (16) 



provided that af/(^ < 1. 

Introducing the shadowing corrections described above into Eq. (|12|), mod- 
ifies it as follows 




We will use this formula in the numerical calculations (see section 3). It should 
be stressed that it is derived in the double log approximation to the DGLAP 
evolution equation. It is instructive to compare this formula with successful 
phenomenology 0, for example, with the Golec-Biernat and Wusthoff approach 
where the ideological basis is shadowing corrections. The Golec-Biernat and 
Wusthoff approach has no impact parameter dependence and thus cannot be 
easily generalized for a nucleus target. Deriving Eq.(20) from AGK cutting 
rules leads to an extra factor exp{—QP'^^) in the second term of Eq.(20) which 
does not appear in the Golec-Biernat and Wusthoff approach. 

For large nuclei the dependence on A of the qqG term is given by the factor 
g-o (x,r^;fex)^ which supprcsscs it especially for large nuclei. 
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3 Numerical calculations 



We now present the results of our numerical calculations with Eq. ([17|). It is 
convenient to carry out the integration over 6^ analytically. We thus obtain 



R 



where 



diff 



I 



dri 



2os(4-)''i J '2 



O"tot 



+ ln(-(fi^^^(x,r^))) 



(20) 



(21) 



/*°*(x,r^;x',rl) = i5;i( V^^(x, r^)) + C + ln(ifi*^^(x, r^)) (22) 



r'^^(x,rx) = 2i?i(fi^^(x,r^) + ^nG(x,r^))-£;i(^^''^(x,r^)) 



MG/ 



+ 21n(-^]G(x,r^) + ^^^'^(x,r^)) 



- ln(fi*^«(x,rx) + l]G(a;,r^)) 



(23) 



where C is the Euler constant, and Ei is the exponential integral function of 
the first kind. 

Our goal in this paper is to study shadowing corrections to DIS. It seems 
reasonable to choose the GRV94 parametrization as input for the DGLAP 
gluon structure function xG^^^^^ . Indeed, GRV94 was derived by fitting ex- 
perimental data at not too small x, so we can hope that the effects of SC which 
are large at small x are not hidden in the specific choice of the initial conditions 
in the parametrization (as is the case in the more advanced global fits to the 
evolution equation |^ HU)- Also, the GRV94 parameterization, with its 
low initial virtuality Qq, is close to the solution of DGLAP in the Leading Log 
approximation. This is compatible with our master equation (0) which has 
been written in the same approximation. 

To avoid double counting we do not utilize more modern parameterizations 
of the structure functions which were determined using data at much lower 
values of x and Q^, and hence are more likely to include some of the screening 
corrections. 
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Before discussing the results of our numerical calculation we would like to 
caution the reader that, we calculate the diffractive dissociation cross section 
integrated over all produced masses, while in the experimental data a mass 
cutoff has been introduced. Calculations, using Eqs. (p0l-|2^), are plotted in 
Figs. 2, 3, 4, 5 and 6. Figs. 2 and 4 show the dependence of R (defined in 
Eq. (21)) on the virtuality . We see that the dependence of R on Q"^ becomes 
steeper as x decreases from x = 10^^ to x = 10"^ over a wide range of 0.5 < 

< 50 GeV^. This implies that we are entering the saturation region, although 
the values that R obtains even for small are far from the unitarity limit 1/2 
(see Eqs.( 0-0)). Obviously, the SC increase with the number of nucleons (A), 
since the parton density scales roughly as A^^^. Figs. 3 and 5 demonstrate the 
increase of the SC with A. 

The same physical behavior is seen in Fig. 6. We see that the energy depen- 
dence of R is far from being flat, contrary to some early expectations 0, which 
were based on the fact that experimentally measured R for a single nucleon is 
flat for mass bins up to 15 GeV[|l|. The main difference of diffractive production 
in DIS on a nucleon target compared to that on a nucleus target is the fact that 
the first depends crucially on the experimental cutoff of the produced mass. At 

~ 1 GeV^ R reaches a value of 0.3 ai W = 220 GeV for DIS on a nucleus as 
heavy as A = 200. Note, that at values of > 10 GeV^ shadowing corrections 
are hardly seen in the HERA kinematical region. 

Finally we draw attention to the Fig. 7 where contributions of qq and qqG 
states to R are shown separately. As seen, the contribution of qqG is almost 
independent of energy, while qq is the main source of energy dependence of R. 

The fact that o"'^*-^-^/crtot is rather small (significantly less than 1/2) for all 
nuclei and it varies substantially with energy, gives rise to a problem of de- 
scribing HERA data within the Kovchegov-McLerran formula. We discuss this 
problem in another publication |T^ . 

We would like to stress that our model completely ignores all corrections to 
the final state (apart from the radiation of gluons). It may turn out that the 
contribution of some of these, e.g. excitation of the nucleus, is rather large and 
can change our results. However it seems unlikely that our conclusions would 
be changed significantly. 



4 Conclusions 

In this letter we studied the influence of shadowing corrections in DIS on nuclei. 
To achieve this goal the generalization of the Kovchegov-McLerran formula ( [T7| ) 
was derived in the Glauber-Mueller approach. This means that we assume that 
a diffractively produced qq state scatters off the nucleons incoherently, and we 
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only keep the Leading Log terms in the formulae for the cross sections. 

Our master equation (|17D makes it possible to calculate the ratio of the 
diffractive dissociation cross section of quark-antiquark pair to the total cross 
section. We take into account the possible contamination of the final state 
caused by the radiation of a single gluon. This contribution is almost indepen- 
dent of energy, and thus we expect that it is qq contribution which saturates 
the ratio R. 

The results of the calculation show that SC start to play a significant role 
in nuclear DIS diffraction in the HERA kinematical region. We thus expect 
that clear experimental signatures of SC will be found in DIS experiments with 
nuclei at HERA and BNL. 
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Figure \.R as a function of the virtuality 
at X = 10~^. From the lower dashed 
curve to the solid one: 
A = 30, 100, 300. 



Figure 5: R as a function of number of 
nucleons A at x = 10~^. Prom the 
upper dashed curve to the solid one: 
Q2 = 1,5, 12, 30, 60 GeV2. 




Figure 6: R as a function of the energy W for different A and Q"^ . From the lower solid 
curve to the upper dashed one: A = 30, 100, 200, 300. 
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Figure 7: Different contributions to R as a function of the energy W for different and 
A. Dashed line: contribution of the qq pair, dotted line: contribution of the qqG state, solid 
line - total contribution. 
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